Disorder-enhanced derealization and local-moment quenching in a disordered 

ant iferromagnet 
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The interplay of disorder and spin-fluctuation effects in a disordered antiferromagnet is studied. 
In the weak-disorder regime (W < U), while the energy gap decreases rapidly with disorder, the 
sublattice magnetization, including quantum corrections, is found to remain essentially unchanged 
in the strong correlation limit. Magnon energies and Neel temperature are enhanced by disorder 
in this limit. A single paradigm of disorder-enhanced derealization qualitatively accounts for all 
these weak disorder effects. Vertex corrections and magnon damping, which appear only at order 
(W/U) 4 , are also studied. With increasing disorder a crossover is found at W ~ U, characterized 
by a rapid decrease in sublattice magnetization due to quenching of local moments, and formation 
of spin vacancies. The latter suggests a spin-dilution behavior which is indeed observed in softened 
magnon modes, lowering of Neel temperature, and enhanced transverse spin fluctuations. 
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I. INTRODUCTION 



The manifestation of quantum antiferromagnetism in 
parent cuprates, JjJ discovered soon after the birth of 
high-T c superconductivity in the doped materials, Q 
has led to intensive efforts to understand the nature of 
this phase both within the Hubbard model, as well as 
its strong-coupling counterpart, the quantum Heisenberg 
model. H Features such as the substantially reduced 
sublattice magnetization (relative to the classical value) 
as deduced from neutron-scattering experiments, jjj the 
substantially enhanced linewidth observed in the two- 
magnon Raman scattering experiments, as well as a 
detailed fitting of the temperature dependence of the spin 
correlation length have confirmed the importance of 
quantum spin fluctuations in these low-dimensional, low- 
spin systems. Antiferromagnetic (AF) spin correlations 
are also present in other strongly correlated systems, no- 
tably the transition- metal oxides such as NiO, V2O3, 
LaV03, NiS2- x Se x , and heavy-fermion compounds such 
as YbP, U2Zni7, UCdn, UCU5 which exhibit AF ordering 
of d and f electrons, respectively, in their ground states. 

Many of the correlated electron systems are intrinsi- 
cally disordered, and the metal-insulator transition ob- 
served in several amorphous materials such as doped 
semiconductors, amorphous Gei_ x Au x and F3i_ x Cu x al- 
loys, and granular alumina have the character of both the 
Mott transition and the Anderson localization transition. 
fl2f The role of strong disorder effects has also been em- 
phasized in the recently-studied transition-metal oxides 
such as LaNii_ x Co x 03 and Na x WC>3. [|l3| The square 
root dip in the electronic density of states near the Fermi 
energy, characteristic of disorder-induced enhancement of 
interaction effects, [[l2| has been found to change to a lin- 
ear form and then to a soft quadratic gap on the insulat- 
ing side. |l4| , |l3| | In some cases, such as in Lai_ x Sr x V03, 
V2O3, the insulator-metal transition is accompanied with 



loss of AF order, [[l3|JlC| l whereas NiS2- x Se x exhibits an 
AF metallic phase HFIl 

It is therefore of interest to study the interplay of quan- 
tum spin fluctuation and disorder effects. Of particular 
interest are questions such as: (i) are transverse spin fluc- 
tuations and quantum corrections to the sublattice mag- 
netization enhanced by disorder? (ii) Does a gapless AF 
state exist, and, if so, at what critical disorder strength 
is AF long-range order (AFLRO) destroyed? (iii) Is the 
AF state destabilized by disorder at finite temperature, 
such that the Neel temperature (for dimensions d > 2) 
is lowered? To answer these questions we will examine 
the influence of diagonal disorder on various properties of 
the Hubbard antiferromagnet, such as sublattice magne- 
tization, quantum spin fluctuations, magnon energy and 
damping, Neel temperature, Hubbard energy gap, and 
electronic density of states. 

Recently spin fluctuation effects were examined 
in impurity-doped antiferromagnets, both within the 
Heisenberg model, [p~5| — p~Q|] as well as the Hubbard model, 
p(j| p2| in order to study magnetic dynamics in cuprate 
antiferromagnets doped with nonmagnetic and magnetic 
impurities such as Zn, Al, Ga, and Fe, Ni, Co respec- 
tively. It was found that a static vacancy, created by 
the replacement of a fermion with a nonmagnetic im- 
purity, for instance, leads to strong magnon scattering. 
It is therefore also of interest to contrast scattering of 
magnons caused by disorder with that by static vacan- 
cies. The third, related, case is that of magnon scattering 
off mobile vacancies, as in hole-doped cuprates, which 
is of course much more efficient at destroying AFLRO; 
spin correlation lengths of order of 1/ ^/x for hole con- 
centration x have been reported from neutron-scattering 
studies. p3fl 

This study therefore complements earlier works on 
the disordered Hubbard model where other aspects have 
been studied, such as the metal-insulator transition, 
p3,E5[ local-moment formation, |26|-|2q] phase diagram 
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etc. A variety of methods have been used earlier, in- 
cluding the scaling theory, field-theoretic approaches, 
p0|- p2| renormalization group (RG ), [|33| , l34| real-space 
RG, [p3p6[ slave-boson formulation, |]37|] dynamical mean 
field theory, []38|-fIo| quantum Monte Carlo studies, ]4l| ] 
and unrestricted Hartree-Fock theory together with ran- 
dom phase approximation, and Onsager-reaction-ficld 
correction to mean-field theory of equivalent spin mod- 
els. [^2|-Q The disordered Hubbard model exhibits ex- 
tremely rich physics and contains the non-interacting 
Anderson localization transition, the purely interact- 
ing magnetic transition, and of course the non-trivial 
fixed point describing the metal-insulator transition in 
the disordered, interacting theory. In addition, vari- 
ous ingredients such as the phenomena of weak localiza- 
tion, disorder-induced enhancement of interaction effects, 
leading to singularities at the Fermi energy, local-moment 
behavior, etc., are also contained. [^5| Generally, the si- 
multaneous presence of interaction and disorder leads to 
a new coupling of the quantum degrees of freedom in 
two-particle quantities that has no counterpart in non- 
interacting, disordered or interacting, pure systems. 

We consider the following Hubbard Hamiltonian with 
random on-site energies, and with a filling of one fcrmion 
per site, so that an AF ground state is obtained, 

H = e ^i<? ~ t /^2 (<*L%> + n - c + U X/^l" - ^' w 

iu (ij)a- i 

The random on-site energies ej are chosen from a uniform 
distribution with —W/2 < €{ < W/2, the distribution 
width W parametrizing the disorder strength. We con- 
sider both the strong correlation limit, with the correla- 
tion term U much larger than the free-particle bandwidth 
B = 2Zt, where Z = 2d is the coordination number, as 
well as the intermediate correlation regime, with U ~ B. 
For concreteness, we consider the square lattice, general- 
ization to three dimension and to other bipartite lattices 
being straightforward. 

We will use several methods/approximations in this 
paper. The broken-symmetry state is obtained in the 
unrestricted Hartree-Fock approximation (UHF), and 
transverse spin fluctuations about this state are studied 
in the Random Phase Approximation (RPA). Disorder 
is treated both perturbatively as well as within a nu- 
merical diagonalization approach on finite lattices. In 
the latter approach the eigenfunctions and eigenvalues 
of the HF Hamiltonian in the fully self-consistent state 
are used to obtain sublattice magnetization, energy gap, 
and the magnon spectrum. The T-matrix approach used 
earlier for impurities, [ p0[ is also employed for compari- 
son. Quantum spin-fluctuation corrections are obtained 
at the one-loop level. ||5| The Neel temperature, energy 
gap, and electronic spectrum are also studied within the 
dynamical mean field theory (DMFT). Q 

The outline of the paper is as follows. Section II deals 
with the reduction of the Hubbard gap due to formation 
of disorder-induced states. Disorder renormalizations of 



the magnon energy, damping, and sublattice magnetiza- 
tion are described in sections III and IV, based on results 
of a perturbative analysis in powers of W/U, discussed 
in the Appendix. A qualitative explanation of the disor- 
der effects is given in Section V in terms of the notion of 
disorder-enhanced derealization. Section VI describes 
the crossover to the strong disorder regime (W > U), 
where the electronic spectrum is gapless, and spin va- 
cancies are formed due to quenching of local moments. 
Magnon softening, enhancement in transverse spin fluc- 
tuation due to spin vacancies, and robustness of the gap- 
less AF state are discussed. Conclusions are given in 
section VII. 



II. DISORDER-INDUCED STATES IN THE GAP 

In this section we examine the formation of disorder- 
induced states within the Hubbard gap, which reduce the 
charge gap, and determine the critical disorder strength 
at which the gap vanishes. We use the T-matrix approach 
and a numerical UHF approach, described in subsections 
A and B respectively. The T-matrix approach is exact 
for a single impurity, and has been used earlier to study 
the formation of defect states within the gap due to a 
single nonmagnetic impurity in the Hubbard AF. p(i| In 
order to use this approach for the disordered AF, with 
random potentials on every site, we make a local approx- 
imation in which we consider a single site, and treat the 
random potential e$ on this site as an impurity potential. 
Comparison with results of the numerical UHF analysis, 
in which disorder is treated exactly, indicates that this 
approximation actually works quite well, particularly in 
the strong correlation limit, where states are strongly 
localized. Spin fluctuation processes will lead to small 
changes in the energy gap of order J in the strong corre- 
lation limit. 



A. T-matrix approach 

Within this approach the location of impurity-induced 
states is obtained from the pole in the T-matrix, T a (oj) — 
e,/(l — ei[g®(w)]u), where [g®]a is the local host Green's 
function. For the pure AF it is given (in the HFA) by, 
= (l/N)YjA?T<rA)/(u 2 -El), for site % in the 
A or B sublattice. Q Here 2A = mU is the Hubbard 
energy gap in the pure AF, = \/ A 2 + e£ is the AF 
band energy, and A is obtained from the self-consistency 
condition (I/AT) Ekl 2 ^) -1 = U' 1 . 

Now, for a given disorder strength W, the inverse po- 
tential (absolute value) has a lower bound l/]e»] > 2/W. 
Therefore poles in the T-matrix are present for |[<7°]m| > 
2/W, so that disorder- induced states are formed within 
the gap, as shown by hatched regions in Fig. 1. If —A 
and A mark the energies (shown by arrows) up to which 
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FIG. 1. Local host Green's function vs - w in the 

pure AF state. Intersections with lines 2/W and — 2/W show 
the extent to which disorder-induced states are formed within 
the Hubbard gap. In all figures energies are in units of the 
hopping parameter t. 

states are formed, then A is obtained from [g®(— A)]u = 
2/W. The remaining Hubbard gap 2A is thus obtained 
from the solution of 



A + A 2 
" W' 



(2) 



A plot of the normalized energy gap 2 A/2 A is shown 
in Fig. 2 as a function of the relative disorder strength 
W/U for U/t— 10, indicating an almost linear reduction 
with disorder strength. 

With increasing disorder strength states are formed 
deeper in the Hubbard gap, and when states have ap- 
proached from both sides in the middle of the Hub- 
bard gap, the energy gap vanishes. The critical disor- 
der strength W c at which the energy gap 2A vanishes 
is therefore given by [g°(0)];i = 2/W c , which yields the 
following equation for the critical disorder strength, 
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A 



2 



(3) 



Considering the strong correlation limit as a special 
case, and keeping terms up to order t 2 /U 2 , the critical 
disorder strength is then given by W c /U = 1 + 8t 2 /U 2 , 



where we used m = 1 — 8i /U and (l/N) X)k e k 
for the square lattice. Thus, with decreasing interac- 
tion strength the ratio W c /U actually increases. This is 
because the kinetic energy becomes relatively important 
with decreasing interaction strength, and the bandwidth 
starts competing with disorder strength. 



B. Numerical UHF analysis 

In order to check the validity of the T-matrix approxi- 
mation, we have also used a numerical UHF analysis. In 
this approach, the HF Hamiltonian on a finite lattice is 
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FIG. 2. Normalized energy gap vs. W/U at T = 0, from 
the T-matrix analysis (line), and the numerical UHF analysis 
for a 10 x 10 lattice (squares). 



numerically diagonalized self-consistently, so that disor- 
der is treated exactly. This approach has been described 
earlier in the context of hole/ impurity doping in the Hub- 
bard antiferromagnet. The energy gap is obtained 
from the energy difference between the lowest energy 
state of the upper Hubbard band and the highest energy 
state of the lower Hubbard band. Configuration averag- 
ing is performed over 100 different realizations of the ran- 
dom on-site potentials on a 10 x 10 lattice. The reduction 
in energy gap with disorder strength using this method 
is also shown in Fig. 2. The numerical analysis shows a 
saturation of the energy difference at W/U ~ 1 due to fi- 
nite system size. Deviations from the T-matrix approach 
are more pronounced at lower interaction strengths where 
the fermion states are more extended. The almost linear 
reduction of energy gap with disorder is also seen at finite 
temperatures, as shown in Fig. 3. The critical disorder 
strength decreases rapidly with increasing temperature. 



III. MAGNON RENORMALIZATION 



Magnons are the low-energy excitations associated 
with transverse spin fluctuations in the broken-symmetry 
state of systems possessing continuous spin-rotational 
symmetry. Therefore they play an important role in 
several macroscopic properties such as the temperature- 
dependence of the order parameter, Neel temperature, 
specific heat, etc. We therefore consider the magnon 
propagator and obtain the disorder-induced renormal- 
izations in magnon energies, transverse spin correlations, 
and the quantum spin-fluctuation correction to sublattice 
magnetization. The magnon propagator with site indices 
i,j is defined in terms of spin-lowering and spin-raising 
operators by GU.+ = (^ G \T[Sr(t)Sf(t')]\^ G ). We take 



the Green's function approach and write the RPA result 
in w-space as 
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FIG. 3. Energy gap vs. W from the numerical UHF anal- 
ysis at different temperatures (shown against plotting sym- 
bols). 



G-+ = 



G° 



(4) 



Here the matrix (in site indices) XcL ' s ^ ne zeroth-order 
antiparallel-spin particle-hole propagator for the disor- 
dered AF, with matrix elements given by [Xdi s ( w )]y = 
i f (du>' /2Tr)\g^(u}')]ij[gi(uj' — uj)]ji, written in terms of the 
one-particle Green's function [ff -( Ct, )]ij> and evaluated in 
the self-consistent, broken-symmetry state. By Xp ure we 
define the corresponding quantity for the pure AF, in 
terms of which the matrix G° = Xpure/ (1 — UXp UIC ) is 
the magnon propagator for the pure AF. Furthermore, 
the disorder self energy E = U 2 8x° is expressed in terms 
of the disorder- induced perturbation 8x = Xdis ~ Xpure- 
Details of the perturbative analysis for 5x° in the strong 
coupling limit are given in the Appendix, and we use the 
result here for the disorder self energy E which has di- 
agonal and nearest-neighbor (NN) terms. For the pure 
AF in the strong coupling limit and in the two-sublatticc 
basis (indices A,B), the propagator in k-space takes the 
form, 



[G°] _1 (k, w) 



1 + w 7k 
7k 1 — u 



(5) 



H,i+S 



1/Z, 



in energy units where 2 J = 8t 2 /U = 1. Here 7k 
: Vos k j + co$k y )/2, so that in real space [G ] -1 has only 
diagonal elements, [G ]^ 1 = 1 ± u> for site i in A and 
B sublattices, and NN matrix elements [G ]^ 1 
where i + 8 refers to NN of i. 

We first consider the configuration-averaged self en- 
ergy (E) = U 2 (8x )- This approximation neglects ver- 
tex corrections which, however, appear only at order 
(W/U) 4 , and are discussed separately in the next sec- 
tion. From the results given in the Appendix for the 
matrix elements of 8x° up to order (W/U) 2 , we obtain 
= -a and (E ii+S ) = -a/Z, in units of U 2 t 2 /A 3 , 
where a = (1/Q)W 2 /U 2 . Substituting the translation- 
ally symmetric (E) in Eq. (4), Fourier transforming 



[G ] -1 — (E), and finally inverting we obtain the follow- 
ing expression for the renormalized magnon propagator 
in k-space and the two-sublattice basis, 



G-+(k,w)=-- 



1 1 



2 V^& 
1 



1 - -ss. 

2J 

~7k 1 
f 



-7k 

2J 



to — ujk + irj lu + d)k — if] 



(6) 



containing both the retarded and advanced parts having 
poles below and above the real-w axis, respectively. Here 
J = J(l + a), and a> k = 2Jyl — 7^ is the disorder- 
renormalized magnon energy. Thus at 0(W/U) 2 the 
form of the magnon propagator remains unchanged, and 
there is only a momentum-independent multiplicative 
renormalization of magnon energies, leading to an effec- 
tive stiffening of the magnon modes by disorder. An up- 
ward shift of the magnon band in the strong-correlation 
and weak-disorder limit has also been observed in a nu- 
merical RPA study in three dimensions. |43| 

This effective enhancement of the magnon energy 
scale can also be viewed as resulting from the enhance- 
ment in the configuration-averaged NN exchange energy 
(t 2 /(U + £j — ej)) which, to second order in W/U, is 
(1 + a)t 2 /U . To the extent that the finite-temperature 
reduction in sublattice magnetization due to thermal ex- 
citation of magnons is suppressed by this enhancement, 
the disordered AF exhibits a slower m(T) vs. T falloff, 
and therefore a higher Neel temperature (for D > 2). For 
strong coupling, the Neel state is therefore stabilized by 
weak disorder, as also reported in the DMFT ffl] and 
the Onsager-reaction- field [H studies. 



Quantum correction to sublattice magnetization 

As only the effective exchange energy scale gets mod- 
ified in Eq. (6), while the form is not changed by 
disorder, |4q | the magnitude of equal-time, same-site 
transverse spin correlations (S^ St) and (S^S^) re- 
main unchanged. These transverse spin correlations are 
obtained by frequency integration of the diagonal ele- 
ments of transverse spin propagators, (S~ (t)S* (£')) = 
— i J(dtu/2iT){G v {uj)]h exp{— iw(t — t')}, where the ap- 
propriate part (retarded or advanced) of G ^ is taken 
depending on whether t' < t or t' > t. Taking t' — ► t~, 
and using the retarded part of G from Eq. (6), we 
obtain 



{Si sf) 



1 1 

2N 



E 



(St Si ) - ^ttE 



k L 



1 

7^1 



+ 1 



(7) 



Here the second result is obtained using the relationship 
Gj^(cj) = Gg^(uj) which follows from spin-sublattice 
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symmetry. The result for total transverse spin fluctua- 
tion (Sfs+ + S+S^), is thus identical to the RPA re- 
sult for the pure Hubbard AF, |Q as well as the spin- 
wave-theory (SWT) result for the Quantum Heisenberg 
Antiferromagnet (QHAF). Therefore up to order 

(W/U) 2 the quantum spin-fluctuation correction to sub- 
lattice magnetization in the strong coupling limit remains 
unchanged from the spin-wave-theory result <5tosf = 
Qi^lEkK 1 - 7k) _1/2 - 1] « 0.39 in two dimensions. 

IIS 



IV. VERTEX CORRECTIONS AND MAGNON 
DAMPING 

In this section we consider vertex corrections which 
were neglected in the previous section due to use of the 
configuration-averaged self energy (E) in Eq. (4). As 
at higher orders this approximation produces terms like 
G°(E)G°(E)G° in Eq. (4), we therefore subtract out this 
term and focus here on the configuration-averaged proper 
self energy at second order, 

r = (EG°E) - (E)G°(E), (8) 

which precisely incorporates the vertex corrections. As 
shown later in this section, this second-order scattering 
process results in magnon damping, which therefore only 
appears at order (W/U) 4 . We note here that the term 
"magnon damping" in this section refers to the decay of 
a momentum mode, and not to the decay into particle- 
hole excitations, which is energetically ruled out in the 
strong-correlation limit where J «U . 

Considering the matrix element , we note that since 
E is only limited to diagonal and NN matrix elements, if 



sites i and j are far apart, then there are no correlations 
between the two E terms, and the difference vanishes. 
Therefore the proper self energy V arises only from local 
correlations in the disorder self energy terms, and has 
diagonal, NN and next-nearest-neighbor (NNN) matrix 
elements. The vertex corrections therefore result in new 
NNN spin couplings in the magnon propagator, besides 
renormalizing the NN couplings. Such longer-range spin 
couplings also arise in the Hubbard AF at intermediate 
and weak couplings. 

For the square lattice (Z = 4), we give below the 
results for matrix elements of T for the case i in A 
sublattice (i G A); results for the other case following 
from symmetry. Also given are the results after sub- 
stitution of the matrix elements of G (k, oj), obtained 
from Eq. (5). We now illustrate the evaluation of T 
for the diagonal matrix element. Expanding the ma- 
trix product, and using the property of the disorder 
self energy that the diagonal element E^ equals the 
sum of the NN elements E^+a, this can be written as 

Tii = [({^i,i+S^i+S' ,i) ~ (^i,i+s) C^i+S' ,i) ) x 

(G% + G° i+S<i + G° i4+S , + G° i+S>i+S ,)}, where summa- 
tion over S and 5' is implied. Configuration averag- 
ing, with cr 2 = {e 2 )/U 2 and a 4 = (ef)/U 4 , for the 
second and fourth moments, yields ((E^ i + s^i+s> i) — 
(E. t .. i+5 )(E l+5Vl )) = Z- 2 [(<T 4 -<T 2 ) + (a 4 + 3a 2 )Sss>}. Sub- 
stituting G°» = (l/7V)E k G°(k, W )exp{zk.(r l - r,-)}, 
and taking the appropriate matrix elements of G°(k,w) 
in the two-sublattice basis, depending on sublattices of 
sites i and j, yields the expression for Fa in Eq. (9). Sim- 
ilarly evaluating the NN and NNN elements, with i + 5, 
i + k, and i + k' standing for the NN, NNN (diagonal), 
and NNN (straight) of i respectively, we obtain, 



Tu = (a* - <J 2 2 )[G% + (G° AB + G° BA ) lk + G° BBl 2 } + Z-\o± + 3a 2 2 )[G° AA + (G° AB + G° BA ) lk + G° BB ] 
= K - a 2 )[(l W )(l - 7 £)/(«£ - u 2 )] + Z-\a 4 + 3af)[2(l j 2 )/(co 2 - c 2 )} 

r M+5 = Z- V - °MG AA + G° BB )(1 + ~f 2 ) + 2(G° AB + G%M + AZ- 2 a 2 [(G AA + G° BB ) + (G% + G BA ) lk ] 
= Z- V " ^2 2 )[2(l ~ Tk)/(^k - + 4Z-V 2 2 [2(1 - 7 2)/(^ _ J>)] 

= Z- 2 (a 4 - (T 2 )2[G° AA cosk x cosk v + (G° AB + G BA )j k + G° BB ] 

= Z- 2 (o 4 - of )2[{2(1 - 7 £) - (1 - w)(l ~ cos k x cosk v )}/(ul - u 2 )] 

T lti+K > = Z~ 2 (o 4 - tT2)[G%72k + (G° AB + G° BA )"/ k + G° BB ] 

= Z- 2 {a 4 a 2 )[{2(l -ft (1 - - 72k)}/(^ - co 2 )} (9) 



Here the summation over momentum k is implied. A straightforward check confirms that the sum of matrix elements 
r,-j + ZTi t i +K + ZTi^ +K > involving sites of the same sublattice (diagonal and NNN) exactly equals the sum Zr^j+a 
involving sites on opposite sublattices. This ensures that the Goldstone mode, which has amplitudes 1 and —1 on the 
two sublattice sites, is preserved, as expected from the continuous spin-rotational symmetry. We also notice that the 
various terms involve fc-sums of the type ^2 k 2 /(c 2 k 2 — to 2 ) from long-wavelength internal magnon modes. Therefore 
the self-energy terms are all non-singular in two dimensions. 
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We now proceed with the magnon renormalization due to this proper self-energy correction T up to order W 4 . To 
this order, it is sufficient to examine the lowest-order correction (q|r|q) to the eigenvalue of the U(l — C^Xpurc) matrix 

for the pure AF. The relevant eigenvalue is 1 — ^Jlo 2 + 7^ in energy units such that 2 J = 1. The magnon amplitudes for 
state |q) are sin 9/2 and — cos 0/2 in the two-sublattice basis, where cos9 = oj / ^Juj 2 + 7^ and sin 6* = 7 q / ^Juj 2 + 7^. 
pl[ For oj = ujq = yjl — 7^ the magnon energy, these amplitudes become yj\ — cj q and — yj\ + w q respectively. Using 
the matrix elements of T from above we obtain for the eigenvalue correction <5A q = (q|r|q), 

(5A q ) = ^ ^[sin 2 (^/2)(r ri + T iti+K Zcosq x cosq y + r^ +K ,Z 72q ) - sin(0/2) cos(0/2)r M+a Z7 q ] 

+ ^ ^ [cos 2 (0/2)^- + I',., . , Z cos 1/, cosg y + T^ j+K ,Zj %tl ) - cos(0/2) sin(0/2)r jJ+5 Z7 q ]. (10) 

In the second term above (for sites j in the B sublattice) the matrix elements of T follow from Eq. (9) with 
u) replaced by — w, in view of Eq. (5). The magnon energy for mode q is then given by the solution of 

1 — yVi 2 + 7 q — <5A q ^ — <5A q \oj) = 0. Here <5A q ^ = — <r(l — 7 q ) is the eigenvalue correction due to the first-order self 

energy (E); its effect on magnon stiffening has been discussed earlier. 

We first consider the magnon renormalization in the long- wavelength (q <C 1), low-energy (to <C 1) limit for 
simplicity. We can drop a-> in the numerators in Eq. (9) for the self energy T, which removes the sublattice dependence, 
and the above eigenvalue correction simplifies in this limit to 

<5A q 2) (o;< 1) = [Tu + T i:i+K Z cos q x cosq y + r i ,i +li ,Z"/2q-T iii+ sZ"/ cl sm9} (11) 
= Z[T lil+5 (l - 7 q sm6») - T M+K (1 - cosq x cosq y ) - I\ i+K /(1 - 7 2q )] 
w aq 2 

As expected the eigenvalue correction goes like q 2 . Here the identity — Z[T i i+ g — T i i+K — T i i+K >] which ensures 
the preservation of the Goldstone mode has been used, and the coefficient a is given by 

a = z [— 2 M+K ' J ' ( ' 

Considering now the case of a general oj, we find that the following terms are present in addition to those given in 
Eq. (11), 

5Xg> - 5\^(u « 1) = ucos8 (1/N) £[(a 4 - a 2 2 )(l - 

k 

— 2Z _1 ((T4 — o\) cos q x cos q y (1 — COS fca, cos fcy)/^ 2 ; — OJ 2 ) 

- Z-\a 4 - ( T 2 2 ) 72q (l - 72k)/(<4 - w 2 )]. (13) 



We now focus on the imaginary part of this second- 
(2) 

order correction 5A q . For this purpose we ex- 
amine the internal momentum sums of the type 
(1/N) ^ k 2uj 2 ./(qj 2 _ — oj 2 ) which appear in the eigenvalue 
correction above. Using the following identity for the 
imaginary part 



^E^^-^E-* k --) 

k k k 

= ttujN^), (14) 

for positive u>, we obtain the imaginary part in terms of 
the magnon density of states N(cu). We consider the two 
limiting cases of low-energy magnon modes w < 1 and 
high-energy modes with u) ~ 1. For long- wavelength, 
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FIG. 4. Spin-dependent electronic densities on the A sub- 
lattice vs. on-site energy in the low-temperature ordered state 
as obtained within the DMFT. The band width is chosen 
B = 8, as in d = 2. 

low-energy modes, N(uS) ~ u>, therefore SMAq 2 '' ~ Dq 2 us 2 , 
where D = ^sa/to 2 in terms of the g 2 coefficient given in 
Eq. (12). We thus find the magnon energy to be given 
by, 

w q = cq[l + iDq 2 ] (15) 

where c is the renormalized magnon velocity. The ratio 
r q /u; q of the magnon damping term to the energy thus 
vanishes like q 2 in the long- wavelength limit, indicating 
weak disorder scattering, caused by the averaging-out of 
the on-site potential disorder at long length scales. 

Long-wavelength magnon modes therefore continue to 
be well-defined excitations even with disorder. However, 
for short- wavelength, high-energy modes with energy 
u> ~ 1, the presence in the imaginary term of the magnon 
density of states, which actually diverges (logarithmically 
in two dimensions) at the upper band edge at energy 2 J, 
indicates that high-energy modes are strongly damped. 
A self-consistent evaluation is therefore required, with 
an imaginary term in the internal magnon mode propa- 
gator. This leads to a self-consistent magnon damping 
of order (W/U) 4 . In as much as the two-magnon Raman 
scattering process probes short-wavelength, high energy 
magnon modes, this magnon damping of order (W/U) 
will be important in an analysis of Raman linewidth in 
disordered antiferromagnets. 

V. DISORDER-ENHANCED DELOCALIZATION 

All the disorder effects in the AF state obtained so far 
can be understood within a single paradigm of disorder- 
enhanced derealization, i.e. an enhancement of the ef- 
fective t/U ratio due to disorder. Thus, the disorder- 
induced enhancement of magnon energy scale and Neel 
temperature in strong coupling, as well as the reduction 
of sublattice magnetization (discussed in the Appendix) , 
can be viewed as arising from this enhancement of the 



FIG. 5. Normalized magnon velocity c(W)/c(0) vs. W/U 
for several interaction strengths, obtained from the low-energy 
magnon modes on a 10 x 10 lattice. 

effective t/U ratio. This derealization effect of disorder, 
contrary to its usual tendency to localize, is character- 
istic of the AF state with its inherent localization due 
to Coulomb barriers, and can be understood as follows. 
A positive on-site energy ti reduces the potential barrier 
for the majority-spin electron, which enhances its proba- 
bility for tunneling through, thereby lowering the on-site 
density. On the other hand, a negative on-site energy 
increases the potential barrier, which has a weaker effect 
if the electron is already localized. Thus the effects of 
positive and negative on-site energies are asymmetrical, 
leading to enhanced derealization on the average. This 
asymmetrical effect is clearly seen in Fig. 4 showing the 
electronic densities for different on-site energies, obtained 
within the DMFT. 

This disorder-enhanced derealization also qualita- 
tively accounts for opposite disorder effects on Tn which 
are observed for strong and weak coupling. While for 
strong coupling, disorder enhances Tn and stabilizes the 
Neel state, the behavior is reversed for weak coupling. 
p0| , [i~4"j Both behaviors can be understood in terms of an 
effective U/t ratio which decreases with disorder. This is 
because in the pure Hubbard antiferromagnet Tn vs. U/t 
actually goes through a maximum, so that with decreas- 
ing U/t, Tn is either enhanced or suppressed depending 
on whether one is on the strong- or weak-coupling side 
of the peak. 

Similar effects are also seen in the magnon velocity, 
c = Wq/g in the limit q — > 0, which also shows a 
peak structure. [|i"6f For this purpose we have consid- 
ered the low-energy magnon modes on a 10 x 10 lat- 
tice. Magnon energies in the RPA are obtained by ex- 
act diagonalization of the x°( w ) matrix, evaluated in 
the self-consistent state using the numerical UHF ap- 
proach. [^6[ The low-energy modes are four-fold degen- 
erate for the pure AF on a square lattice, and therefore 
the normalized magnon velocity is obtained by averaging 
Lo n (W)/u} n (0) over the four lowest energy modes. The 
normalized magnon velocity, configuration-averaged over 
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20 different configurations, is shown in Fig. 5, for sev- 
eral values of U/t. For strong coupling, the magnon ve- 
locity increases with disorder, almost quadratically, in 
agreement with the perturbative result from Eq. (6). 
The behavior is reversed for weak coupling, and the 
magnon velocity decreases with disorder strength. As 
the temperature-dependence of the sublattice magneti- 
zation in the low-temperature limit depends only on the 
energy scale of low-energy magnon modes, this also im- 
plies slower (faster) temperature-falloff of m in the strong 
(weak) coupling limit. 



VI. STRONG DISORDER (W > U) AND 
CROSSOVER BEHAVIOR 



We have seen earlier that with increasing disorder 
strength, the energy gap decreases almost linearly with 
W, and eventually closes when W ~ U. Therefore for 
W > U, assuming a single-occupancy constraint, the 
two Hubbard bands would overlap, indicating that sin- 
gle occupancy for all sites is no longer energetically fa- 
vorable. Electrons from the highest-energy sites (with 
€. L > £7/2) are transferred to the lowest-energy sites 
(with < —U/2), making them doubly occupied and 
shifting up the energy by £7. The electronic states as- 
sociated with these essentially empty and doubly occu- 
pied sites are therefore located near the Fermi energy at 
£7/2. As these unoccupied and doubly occupied sites are 
nonmagnetic, the sublattice magnetization starts falling 
relatively rapidly with disorder strength. In the limit 
t <C £7, W the bands are nearly flat with bandwidth ~ W, 
and the overlap region is ~ (W — £7), so that a simple 
estimate for the fraction of these nonmagnetic sites yields 



W-U 

w 



(16) 



indicating an almost linear decrease in sublattice mag- 
netization for W ~ U . Thus a crossover takes place at 
W ~ £7 from the essentially flat sublattice magnetization 
to an almost linear falloff with disorder. This is clearly 
seen in Fig. 6 where the configuration-averaged sublat- 
tice magnetization obtained within the UHF approxima- 
tion is plotted against the disorder strength W. 

However, a more significant consequence of the forma- 
tion of these nonmagnetic sites is that they essentially act 
like spin vacancies in the antiferromagnet, which leads to 
spin-dilution behavior, as discussed in the following sub- 
section. The above picture suggests that for W > £7 the 
system can be viewed as a composite of a disordered AF 
with W = £7, and a spin-diluted system with a concentra- 
tion x = (W — U)/W of spin vacancies. Strong magnon 
scattering off static vacancies, leading to substantial soft- 
ening of low-energy, long-wavelength modes and magnon 
damping has been obtained earlier. [ ^6|Jl^j2l|] The ra- 
tio of magnon damping term to its energy now goes like 
(W — U)/W, to be contrasted with the small damping 



ratio (W/U) q 2 in the weak-disorder regime, obtained 
earlier in section IV. Vacancy-induced enhancement in 
transverse spin fluctuation has not been studied, and we 
discuss this in the following subsection. 



A. Enhanced fluctuations due to spin vacancies 

We present a simple estimate of the enhancement of 
transverse spin fluctuations due to spin vacancies. We 
consider the following Hubbard Hamiltonian on a square 
lattice with binary-distributed, random NN hopping, 



H 



E 

<ij>0 



tiA&ia&jv +h.C.) + £7 y^nii^ix 



(17) 



where the hopping term — if sites i or j are vacancy 
sites, and Uj — t otherwise. Thus for a vacancy on site 
i, all hopping terms U^+g connecting i to its NN sites 
i + S are set to zero. The vacancy site is thus completely 
decoupled from the system. Half-filling is retained by 
having one fermion per remaining site. We consider the 
U/t — ► oo limit, where the local moments are fully satu- 
rated, and the model maps to the localized-spin Heisen- 
berg model. In this limit the vacancy problem becomes 
identical to the spin vacancy problem in the QHAF, for 
which magnon renormalization was studied earlier. p6[ 

The structure of the x ( w ) matrix in the host AF, and 
the modification introduced by spin vacancies has been 
considered earlier in the context of static impurities. 0] 
Since the vacancy spin is completely decoupled from the 
system, the magnitude of the diagonal matrix element 
[x ]u on the vacancy site i is irrelevant. To minimize the 
perturbation, we treat the vacancy site as occupied with 
an isolated spin. For a vacancy on site i, a perturba- 
tion is induced in the neighborhood due to the absence 
of hopping between vacancy site and NN sites i + 6. In 
terms of the notation used in section III, the following 
self energy correction E = U 2 8x° is obtained 



Ejj — 1 



1/4 



(18) 



where again 2 J has been set to 1. Thus, a diagonal con- 
tribution Eji arises if a vacancy exists either on site i, 
or on any of the four NN sites i + 6. For a finite va- 
cancy concentration x, the probability that a vacancy 
exists on a site is x. Therefore, configuration-averaging 
yields (£«) = 2x, and (E, M+(5 ) = 2ar/4. From Eq. (4) 
we obtain the same expression for the magnon propa- 
gator as in Eq. (6), except that now J = J(l — 2x), 
so that d>k = 2J(1 — 2x) \/l — 7^ yields a softening of 
the magnon mode, reflecting the spin-dilution behavior. 
pel This magnon softening is in contrast to the stiffening 
in the weak-disorder regime, obtained earlier in section 
III. Enhanced thermal excitation of magnons due to this 
softening will result in a faster temperature falloff of the 
sublattice magnetization, and hence to a lowering of the 
Neel temperature. 
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FIG. 6. Sublattice magnetization m vs. W, at the HF level FIG. 7. Neel temperature T N vs. W obtained within the 

for a 10 x 10 lattice (dashed line), and including spin fluctu- DMFT 
ations (solid line). 



Since the form of the magnon propagator is not 
changed, as already discussed in section III, the trans- 
verse spin correlations (S~ S^) and (S^~S~), as well as 
the lattice-averaged spin-fluctuation correction to sublat- 
tice magnetization remain unchanged. Therefore, from 
the SWT result given after Eq. (7), Sm SF w 0.39 in two 
dimensions. However, as the number of spins is reduced 
to 1 — x, the quantum correction per spin is enhanced. 

In a recent numerical, finite-size study of the spin- 
vacancy problem, focussing on the exact eval- 
uation of transverse spin fluctuations in the RPA, 
the x-dependence of the quantum correction per spin 
<5msF/(l — x) was found to be best described by the ex- 
pression 0.39 + 0.42a; + 6.5a; 3 . This results in a nearly 
vanishing O(x) term in <5mgF, in agreement with the 
perturbative result obtained above. The lattice-averaged 
sublattice magnetization for the disordered AF, obtained 
by accounting for the quantum spin-fluctuation correc- 
tion using m — hihf — Sm^F, with SmsF — (1 ~~ 
a;)(0.39 + 0.42a;-|-6.5a; 3 ), where the vacancy concentration 
x = (W - U)/W for W > U, is compared with the HF 
result in Fig. 6. The crossover at W ~ U is again clearly 
seen from the rapid decrease in sublattice magnetization. 



B. Gapless antiferromagnetic state 



diagram of which has been recently studied within the 
UHF approach. [Q A region of metallic AF state in d = 3 
was also recently identified in the t — i'-Hubbard model 
without disorder; here it is the NNN hopping amplitude 
t' that leads to frustration of AF order. j5^ ] 

The robustness of the gapless AF state is also seen from 
a DMFT calculation of the Neel temperature Tn, for a 
rectangular distribution of e^, as obtained from an ex- 
trapolation to zero of the inverse AF susceptibility. Both 
exact enumeration and Monte Carlo evaluation were em- 
ployed. Details of the application of DMFT to the dis- 
ordered Hubbard model have been described earlier, pCj ] 
where binary alloy and semicircular distributions were 
studied. The variation of Tn with W is shown in Fig. 7 
for U/B = 1. This corresponds to U/t = 8 in the two- 
dimensional case. Appreciable spin ordering is evident 
from the fairly high Tn even at W ~ U where the energy 
gap vanishes. Within the DMFT, this closing of the en- 
ergy gap with disorder, in the low-temperature ordered 
state is seen in Fig. 8 from the single-particle density 
of states N(ui), which was obtained by analytically con- 
tinuing the imaginary-time Green's function using the 
maximum entropy procedure. [ p3[ 



The above result indicates that substantial AF order- 
ing remains even for W ~ U, where the energy gap 
closes. Since AF order persists for W > U the gapless 
AF state is apparently quite robust even in two dimen- 
sions. This feature was also observed in the QMC study 
in d = 2, |ll|] where substantial AF correlations were seen 
for W = U = 4t, while the compressibility indicated an 
absence of the charge gap. This leads to the possibil- 
ity, in three dimensions, of a metallic AF state, if states 
at Fermi energy are not localized. Gapless AF states, 
both metallic and insulating, were also obtained for the 
three-dimensional disordered Hubbard model, the phase 



As expected for weak and intermediate couplings, Tn is 
seen to decrease with disorder strength, in contrast to the 
strong-coupling result of an enhancement in the magnon 
energy scale, and hence in Tn- Also, the critical disorder 
strength where Tn vanishes is seen to be nearly 1.5 U, 
which is intermediate between the critical values of nearly 
2 U for the milder semicircular distribution and nearly U 
for the much more severe binary distribution. |io| In fact, 
for the binary alloy case, the zero-temperature transition 
is expected to occur at W « U, when all sites abruptly 
become either unoccupied or doubly occupied and the 
sublattice magnetization vanishes. 
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FIG. 8. Electronic density of states N(u) in the 
low-temperature ordered state for several disorder values, 
as obtained within the DMFT. The charge gap closes near 
W = U. 



VII. CONCLUSION AND DISCUSSION 

Two fundamentally different mechanisms — disorder- 
enhanced derealization and local-moment quenching — 
were identified to control the magnetic behavior of an AF 
in the regimes of weak and strong disorder, respectively. 
In the weak disorder regime (W < U) disorder effects 
on sublattice magnetization, magnon-mode energies, and 
Neel temperature can be qualitatively understood within 
the disorder-enhanced delocalization effect. The stabi- 
lization of the Neel state by disorder in the strong corre- 
lation limit, reflected in an enhancement of the magnon 
energy scale and the Neel temperature, is a striking con- 
sequence. In this regime the AF state is remarkably ro- 
bust against disorder, particularly in the strong correla- 
tion limit. Low-energy, long- wavelength magnon modes 
are weakly damped, and continue to be well-defined ex- 
citations. 

With increasing disorder strength there is a crossover 
at W ~ U, characterized by a rapid decrease in sublat- 
tice magnetization and quenching of local moments due 
to formation of nonmagnetic sites. Driven by band over- 
lap, this local-moment quenching can be viewed as the re- 
verse of local-moment formation in the disordered metal- 
lic state with increasing U /W. These nonmagnetic sites 
act like spin vacancies in the antiferromagnet, leading to 
characteristic spin-dilution behavior. Vacancy-induced 
magnon scattering results in enhanced transverse spin 
fluctuations, softened magnon modes, lowering of Neel 
temperature, and strong magnon damping. 

As discussed in section VI, in the gapless AF state 
for W > U, the electronic states near the Fermi energy 



are associated with (and localized around) the essentially 
empty and doubly occupied sites. These states are non- 
magnetic in that they do not contribute significantly to 
local moments. Therefore it is interesting to note that 
the low-frequency conductivity behavior involves these 
nonmagnetic states near the Fermi energy. Further in- 
vestigation of the gapless AF state, focussing on the low- 
frequency conductivity behavior, magnon damping due 
to decay into particle-hole excitations, and an estima- 
tion of the critical disorder strength where AFLRO is 
destroyed, is presently in progress. p4j | 
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APPENDIX 

PERTURB ATIVE ANALYSIS IN W/U: STRONG 
CORRELATION LIMIT 

In this appendix we treat the disorder term V = 
Y] il7 eihi a in the Hamiltonian as a perturbation and ob- 
tain disorder corrections in the broken-symmetry state 
up to second order in W/U. With g Q a representing the 
one-electron HF Green's function in the pure AF, the cor- 
responding Green's function for the disordered AF, g a — 
g®/[l — Vg®}, then yields disorder corrections to elec- 
tronic densities and sublattice magnetization. Similarly, 
the magnon renormalization within RPA is obtained from 
the correction dx° = Xdi s — Xpurc m the zeroth-order 
antiparallel-spin particle-hole propagator defined below 
Eq. (4). 

For analytic convenience we consider the strong- 
correlation limit and retain terms of order t 2 /U 2 only, 
at which level the pure AF is equivalent to the NN 
Heisenberg model. Up to this order the Green's function 
[g®]ij in site-basis has only diagonal, nearest-neighbor 
(NN) and next-nearest-neighbor (NNN) matrix elements. 
Only the former two are actually required in the strong- 
coupling analysis, and are given below for spin up and 
site i in the A sublattice. Expressions for other spin 
and sublattice cases follow from the spin-sublattice sym- 
metry. These Green's functions are easily obtained by 
starting with the atomic limit, where the bands are at 
energies —A and A, with 2 A = U, and then obtain- 
ing corrections up to second order in the hopping term 
t. The diagonal and NN matrix elements of the time- 
ordered Green's function matrix [g®], containing both the 
advanced (lower band) and retarded (upper band) parts 
are, 
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[9? 



l-i 2 /A 2 



i 2 /A 2 



[fl , t]*,i+<5 



J " ~w-(-A)-tt7 
-t/2A 



w — A + z?7 
i/2A 



cj — (—A) — i?7 w — A + irj 



(19) 



We note here that the corrections to one-particle 
Green's functions involve renormalizations of (i) density 
(wavefunction) and (ii) energy. Diagrams in which the 
potential scattering is interband (involving different Hub- 
bard bands) result in transfer of spectral weight across 
the Hubbard bands, and therefore yield density changes. 
However, the intraband processes (involving all propaga- 
tors from the same band) represent energy renormaliza- 
tion due to disorder potential, and do not contribute to 
any density change. In the particle- hole propagator x°, 
while both energy and density renormalizations need to 
be considered, diagrams involving only energy renormal- 
ization do not contribute as both the particle and hole 
energies are shifted identically by the disorder potential. 



1. First order 

The first-order correction g®Vg® yields the following 
interband contribution to the local one-particle Green's 
function for spin up, 



[^9^ \ii I interband 2 

-t/2A 



1 



t 2 /A 2 



uj — (—A) — ir] ui — A + ir] 
t/2A 



f uj - (-A) - ir] 



UJ — A + IT] 



(20) 



where i + 6 refers to the NN sites of i. The correction 
to density results from the spectral weight transferred to 
the upper band, 



5n$ = 



Ml.. - 



A 2 



U 



Thus for positive q (and 6i+$ = 0, for sake of argu- 
ment) the electron on site i is more delocalized as its 
energy is pushed up, leading to the above decrease in 
density. A simple way to see this is in terms of the es- 
cape probability for the spin-up electron from site i to its 
nearest neighbors i + 5 due to the virtual hopping pro- 
cess. The net probability of escape changes from t 2 /U 2 
to (1/4) J2s t 2 I (U + e i+s ~ e i) 2 ) which to first order in e 
precisely yields the above result. 

From the particle-hole symmetry with spin flip in the 
Hubbard antiferromagnet, fiq| the spin-down hole expe- 
riences the same effect, except that the potentials are 
reversed for the hole. Therefore, there is an increase of 
identical magnitude in the spin-down hole density, im- 
plying an identical decrease in the spin-down particle 
density. This leads to a vanishing change in the local 
magnetization tohf, and the particle- hole propagator \ ■ 



This cancellation would be absent if the random potential 
was spin-dependent, and is connected to the time-reversal 
symmetry. Thus, 



5m 



(i) 

HF 



0; [5 X °] (1) =0 (22) 



2. Second order 

Second-order changes in densities, Sn ia , are similarly 
obtained from the second order correction g®Vg®Vg®. As 
this term is invariant under the transformation V — > — V, 
identical changes are obtained for spin-up electron and 
spin-down hole. The changes in spin-up and spin-down 
electronic densities are therefore equal and opposite, and 
hence the change in local magnetization adds up to, 



8m 



(2) _ af_l V ( e * - e i+s) 2 
" A 2 4V U 2 

6 



(23) 



This result also follows from the second-order correc- 
tion to the net escape probability, discussed earlier below 
Eq. (22). Configuration averaging, with (e 2 ) = W 2 /12, 
therefore yields the following quadratic decrease of sub- 
lattice magnetization with disorder, which is suppressed 
by the small factor t 2 /U 2 . As discussed in section 
III, quantum spin-fluctuation effects do not substantially 
modify this result, which qualitatively agrees with the 
quantum Monte Carlo studies of the disordered Hubbard 
model. Bill 



u 2 u 2 ' 



(24) 



Turning now to the second-order corrections to x > we 
find that processes containing energy renormalization in 
the fermionic propagators (either solely, or along with 
density renormalization) cancel, so that the net result to 



(2i) order {W/Uf 



0l(2) 



MS* 



A 3 4 



6 

A 3 4 



u 2 

- e i+s) 2 
U 2 



(25) 



We note that the sum of all NN matrix elements is pre- 
cisely the diagonal matrix element. An immediate con- 
sequence of this correlation is that the Goldstone mode 
is preserved, as expected from spin-rotational symmetry, 
and that generally the effective scattering of low-energy, 
long-wavelength magnon modes is weak. This disorder- 
induced perturbation 5x° directly yields the magnon self 
energy, and disorder renormalization effects on magnon 
properties are discussed in sections III and IV. 
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